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SENSIBLE  AND  LATENT  HEAT  FLUXES  AND 
HUMIDITY  PROFILES  FOLLOWING  A  STEP 
CHANGE  IN  SURFACE  MOISTURE 


Edgar  L.  Andreas 

INTRODUCTION 

The  velocity,  temperature,  and  humidity  profiles 
of  Rider,  Philip,  and  Bradley  (1963),  hereafter  re¬ 
ferred  to  as  RPB,  are  an  ideal  set  for  testing  models 
of  air  flow  encountering  an  abrupt  change  in  surface 
conditions  and  have  been  used  several  times  for  this 
purpose  (Townsend  1965,  Taylor  1971,  Rao  et  ai. 
1974).  Briefly,  RPB  studied  the  modification  in  a 
flow  blowing  from  a  hot,  smooth  tarmac  airport  run¬ 
way  onto  a  well-watered,  grassy  field  (Fig.  1 )  at 
Canberra,  Australia.  Their  data  set  includes  tempera¬ 
ture  and  humidity  at  five  heights,  between  5  and  1 50 
cm,  upwind  and  at  several  fetches  downwind  of  the 
surface  change,  a  velocity  profile  that  is  the  average 
of  the  upwind  profile  and  the  downwind  profile  16 
m  from  the  surface  change,  and  measurements  of 
upwind  and  downwind  net  all-wave  radiation,  ground 
heat  flux,  and  radiative  surface  temperature. 

Their  data  set  is  so  complete  that  RPB  could  re¬ 
port  on  only  a  fraction  of  the  possible  computations 
it  allows.  Sibbons  (Rider  et  al.  1 965)  pointed  out  a 
few  of  their  omissions.  The  most  obvious  are  calcu¬ 
lations  of  the  upwind  and  downwind  sensible  and 
latent  heat  fluxes.  RPB  did  discuss  the  average  down¬ 
wind  latent  heat  flux  but  did  not  report  fluxes  for 
individual  runs.  Thus,  when  Rao  et  al.  (1974)  used 
the  data  set  in  their  numerical  model,  they  based 
surface  flux  estimates  on  the  net  radiation  and  ground 


flux  values  that  RPB  did  report.  In  testing  his  model, 
however,  Taylor  (1971)  calculated  the  upwind  sen¬ 
sible  heat  flux  for  two  runs  from  the  given  profiles 
and  found  these  values  smaller  than  the  sum  of  the 
net  radiation  and  ground  heat  flux.  Another  look  at 
the  RPB  data  set,  with  emphasis  on  computing  the 
surface  heat  fluxes  from  the  reported  profiles,  would 
therefore  be  valuable. 

I  also  have  two  personal  reasons  for  looking  at  the 
data.  Andreas  et  al.  (1979)  measured  the  sensible 
heat  flux  over  arctic  leads  in  winter  and  found  that 
this  flux  obeyed 

Ni  -  0.14  Rx0-72,  (1) 

where  Ni  is  the  sensible  heat  flux  Nusselt  number  and 
Rx  is  the  fetch  Reynolds  number.  Their  attempt  to 
measure  latent  heat  flux  was  unsuccessful,  however; 
they  simply  had  to  assume  that  temperature  and 
moisture  were  transferred  similarly-that  the  latent 
heat  flux  Nusselt  number  also  obeys  eq  1 .  In  the 
flow  RPB  studied  there  was  a  large  step  increase  in 
latent  heat  flux  downwind  of  a  uniform  surface  with 
negligible  flux-just  what  occurs  over  leads.  I  there¬ 
fore  hope  that  the  RPB  data  will  help  us  better  under¬ 
stand  heat  transfer  from  leads. 

The  high  quality  RPB  humidity  profiles  also  in¬ 
terest  me.  Andreas  et  al.  (1981)  reported  that  over 
arctic  leads  the  condensate  concentration  profile, 

C(z),  obeys  the  relation 


Figure  1.  Internal  boundary  layer  over  the  downwind  grass. 


l-CW/c.|*  =  -0.34lnU/A),  (2) 

where  c«  is  a  concentration  scale  formed  from  the 
condensate  flux,  and  h  is  an  estimator  of  internal 
boundary  layer  height  that  depends  on  fetch.  Since 
this  was  the  first  report  of  such  a  profile  form,  we 
seek  confirmation.  Andreas  etal.  (1981)  inferred 
that  the  form  of  eq  2  was  distinctly  different  from 
that  of  the  temperature  profiles  over  leads  because 
condensate  was  transferred  passively  while  tempera- 
was  not.  Because  humidity  was  similarly  passive  in 
the  regime  RPB  studied,  I  look  for  evidence  in  their 
downwind  humidity  profiles  to  support  eq  2. 


MATHEMATICAL 

Upwind:  flux-gradient  method 

The  surface  energy  budget  for  steady-state  condi¬ 
tions  is 

fln-G  =  //s  +  WL.  (3) 

Here  I  have  followed  Munn’s  (1966)  sign  conven¬ 
tion:  a  gain  in  surface  energy  from  net  radiation,  Rn, 
is  positive;  a  ground  flux,  G,  into  the  ground  is  posi¬ 
tive;  the  turbulent  sensible  (Ws)  and  latent  (WL)  heat 
fluxes  are  both  positive  when  upward.  RPB  reported 
Rn  both  over  the  tarmac  and  1 6  m  downwind  over 
the  wet  grass.  They  measured  G  at  these  same  loca¬ 
tions  for  25  of  their  43  experimental  runs.  I  will 
calculate  H%  and  H L  from  their  velocity,  temperature, 
and  humidity  profiles. 

The  calculations  are  complicated  because  RPB  re¬ 
ported  a  velocity  profile  that  is  the  average  of  upwind 
and  downwind  profiles.  From  their  Figure  5  and 
Figure  1  in  Rider  et  al.  (1965),  the  downwind  veloc¬ 
ities,  especially  in  the  lowest  meter,  clearly  were 
lower  than  upwind  velocities.  Using  their  average 
profile  would,  therefore,  bias  both  upwind  and  down¬ 
wind  flux  calculations:  I  need  individual  upwind  and 
downwind  profiles.  Because  the  two  figures  men¬ 
tioned  show  that  upwind  and  downwind  profiles  were 
virtually  identical  for  heights  above  1  m,  the  individ¬ 
ual  profiles  can  be  retrieved  from  the  average  velocity 
profiles  RPB  tabulated  by  the  flux-gradient  method 
(Businger  et  al.  1971 ).  I  simply  assume  that  over  the 
uniform  tarmac  both  upwind  velocity  and  tempera¬ 
ture  profiles  obeyed  Monin-Obukhov  similarity  and 
fit  them  with  empirical  similarity  functions. 

There  admittedly  is  still  no  universal  agreement 
over  the  forms  of  these  similarity  functions  (Yaglom 
1977),  but  Wieringa's  (1980)  recent  analysis  of  the 
effects  of  flow  distortion  on  the  Kansas  results 


(Businger  et  al.  1971  )-the  most  discordant  of  the 
published  similarity  functions-has  resolved  the  main 
points  of  controversy.  For  the  unstable  flow  upwind, 
I  use  the  integrated  similarity  functions  from  Paulson 
(1970): 

U{z)  [ln(z/z0)  -  'bmklL)\  (4) 

and 

T{z)=  r01  +  0.[ln(z/zh)-^h(z/Z)|,  (5) 

where 


vl'm  =  2  ln[(1+y)/2]  +  ln[(1+y2)/2] 


-  2  arctan  y  +  n/2, 

(6) 

^h  =  2ln[(1+y2)/2], 

(7) 

y=  (1-16 (zlL)\v‘  . 

(8) 

The  Obukhov  length  corrected  for  moisture  flux 
(Zilitinkevich  1966,  Busch  1973),  a  measure  of  sur¬ 
face-layer  stability,  is 


u*3 

9  k  wOj 


x 


o.6i  r,7 

i  + - — 

L  1  +  0.61  Q- 


wdj 


(9) 


I  fitted  eqs  4  and  5  with  the  RPB  upwind  temperature 
data  at  z  =  5,  1 1 .5,  27.5,  64, 1 50  cm  and  with  the 
velocity  data  at  z  =  1 00, 1  ;>0,  200  cm— levels  that  I 
believe  reflect  the  upwind  velocity-and  at  z0,  where 
U(z0)  =  0. 

The  flux-gradient  method  is  an  iterative  procedure. 

I  initially  assume  neutral  stability —that  is,  z/L  =  0, 
then  =  i)/h  =  0— and  make  a  least-squares  fit  of 
the  data  compatible  with  the  models,  eqs  4  and  5. 

The  slopes  of  the  least-squares  lines  yield  estimates 
of  ut  and  0,  from  which  I  can  estimate  z/L  and,  in 
turn,  and  i//h.  Using  these  and  iph  values  in 
eqs  4  and  5,  I  again  pass  least-squares  lines  through 
the  data.  More  accurate  values  of  u,  and  0,  result. 

The  process  continues  until  the  relative  sum  of  squared 
profile  deviations  changes  by  less  than  0.5%  between 
consecutive  iterations.  This  usually  occurs  after  no 
more  than  five  iterations. 

RPB  gave  0.002  cm  as  the  upwind  value  of  z0, 
Taylor  (1 971 )  questioned  the  smallness  of  this  number, 
however;  with  z0  =  0.002  cm,  the  roughness  Reynolds 
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number,  R,  ~  u,z0jv,  is  0.2  -  0.3,  barely  above  the 
upper  limit  for  oerodynamically  smooth  flow,  0.1 3, 
and  well  below  the  lower  limit  for  oerodynamically 
rough  flow,  2.5  (Businger  1973).  In  using  z0  =  0.002 
cm  I  find  that  the  flux-gradient  iteration  does  not 
accurately  reproduce  the  upwind  and  downwind 
profiles  shown  in  Figure  5  of  RPB  or  Figure  1  of  Rider 
et  al.  (1965).  The  value  z0  =  0.01  cm,  a  roughness 
length  consistent  with  ones  reported  for  similar  sur¬ 
faces,  yields  upwind  and  downwind  velocity  pro¬ 
files  that  agree  well  with  the  available  profiles  (i.e., 
Figure  5  of  RPB  and  Figure  1  in  Rider  et  al.  1965). 

In  eq  5  ^01  is  assumed  to  be  the  surface  tempera¬ 
ture  of  the  tarmac.  After  calculating  6,  we  can  esti¬ 
mate  T0t  by  assigning  a  value  to  irh,  the  roughness 
length  for  temperature.  I  find  that  with  zh  =  0.2 
cm,  Tg,  values  computed  with  eq  5  are  usually  with¬ 
in  2°C  of  the  radiative  temperatures  measured  by 
RPB  and  also  agree  well  with  Sibbon’s  (Rider  et  al. 

1965)  calculations  of  Tot  for  runs  6  and  7. 

These  estimates  of  r01  permit  evaluation  of  the 

downwind  surface  temperatures  T02  from  the  net 
radiation  values  tabulated  by  RPB.  I  simply  decom¬ 
pose  the  upwind  and  downwind  values  of  net  radia¬ 
tion  into  shortwave  and  longwave  components  (Munn 

1966) : 

K„,  =(1  -a])Rs*RL-e1or0J4,  (10) 

/?n2  =  (1  -a2)/?s  +  «L-e2ar024.  (11) 

In  these  Rs  is  the  incident  shortwave  radiation;  Oj, 
the  surface  albedo  (/'  =  1,2);  RL,  the  incident  long¬ 
wave  radiation;  and  CjOT^j4,  the  emitted  longwave 
radiation.  RPB  stated  that  within  the  accuracy  of 
their  measurements  a,  =  a2  •  Hence,  eqs  1 0  and  1 1 
show  that 

T02  =|(^nl  ~Rn2^°€2 

♦(‘iWoi4!*-  02) 

With  realistic  values  for  the  emissivity,  =  0.95 
and  e2  =  0.99  (Kondratyev  1969,  Siegel  and  Howell 
1971 ,  Razurnovskii  and  Popov  1974),  eq  12  yields 
values  of  T02  8°C  lower,  on  the  average,  than  the 
radiative  temperatures  RPB  listed.  My  T02  values 
for  Runs  6  and  7  do,  however,  agree  very  well  with 
Sibbon’s  estimates  (Rider  et  al.  1965).  The  down¬ 
wind  temperature  and  humidity  profiles  imply  that 
my  values  of  f02  are  more  reasonable  than  those 
RPB  tabulated.  Their  downwind  radiative  tempera¬ 
tures  are  usually  higher  than  the  temperature  at  the 
lowest  downwind  profile  level.  This  not  only  neces¬ 


sitates  a  complex  shape  for  the  temperature  profile 
near  the  surface  but  also  implies  a  positive  sensible 
heat  flux  over  the  grass,  a  result  incompatible  with 
the  shape  of  the  downwind  temperature  profile. 

When  compared  with  several  RPB  profile  sets,  the 
numerical  models  of  Taylor  (1971)  and  Rao  et  al. 
also  predicted  downwind  surface  temperatures  con¬ 
sistently  lower  than  those  RPB  listed  in  their  Table  2. 

The  downwind  humidity  profiles  confirm  the 
accuracy  of  my  downwind  surface  temperature 
values.  In  the  heat  flux  calculations  that  follow  I 
use  for  the  downwind  surface  humidity  <?02  the 
value  resulting  from  an  extrapolation  according  to 
eq  37  through  the  lowest  two  profile  levels  to/02  = 
0.14  cm.  Andreas  et  al.  (1981)  used  this  same  ex¬ 
trapolation  to  calculate  surface  condensate  concen¬ 
trations  over  leads.  The  resulting  values  are  always 
near  the  humidities  computed  assuming  saturation 
with  a  surface  at  temperature  T02.  If,  on  the  other 
hand,  the  RPB  surface  temperatures  were  used,  the 
surface  relative  humidity  would  generally  be  about 
70%,  contrary  to  the  assumption  by  RPB  and  Sibbons 
(Rider  et  al.  1965)  that  the  air  is  saturated  over  the 
wet  grass. 

Downwind:  integral  method 

An  internal  boundary  layer  (IBL)  developed  when 
air  that  was  essentially  in  equilibrium  with  the  hot, 
dry  tarmac  runway  encountered  the  cool,  moist 
grass  (Fig.  1 ).  Air  within  the  IBL  was  modified  by 
the  fluxes  from  the  new  surface,  while  the  air  above 
was  unmodified.  The  existence  of  an  IBL  is  the 
basis  for  my  calculation  of  the  downwind  fluxes. 

Consider  a  control  volume  over  the  downwind 
grass  as  shown  in  Figure  1 .  The  top  of  this  volume 
is  atz  =  z  ,  where  z  is  greater  than  the  IBL  height 
5fx)  for  all  x.  Andreas  etal.  (1979, 1981)  have  used 
the  fact  that  the  heat  and  mass  content  of  this  volume 
is  constant  in  steady-state  conditions  to  derive  rela¬ 
tions  between  surface  fluxes  and  profiles  measured 
at  the  upwind  and  downwind  ends  of  the  control 
volume.  For  completeness,  I  will  go  through  an 
abbreviated  derivation  here;  to  fill  in  gaps,  the  reader 
can  refer  to  Andreas  et  al.  (1979, 1981 ). 

The  heat  budget  of  the  control  volume  is 

pf  U^znc^z)^  LvQ.}(z))  dz 

J  o 


-pf  M/(x,z’)[cp7'i(/)  +  Lv(?j(z’))  dx 
J0 


3 


Consequently,  the  sensible  and  latent  heat  budgets 
of  the  volume  are  separable: 


rx  —  — 

-PI  (cpw0j  +  Lvwqt)dx 

J0 

-fif  U(X,2)[cpnx,2) 

J0 

+  L,Q(X,z)\  dz  =  0.  (13) 

Here  U,  W,  T,  and  Q  are  average  values  of  longitudinal 
and  vertical  velocity,  temperature,  and  specific  hu¬ 
midity;  w,  8,  and  q  are  the  turbulent  fluctuations  of 
these;  p  and  cp  are  the  density  and  specific  heat  of 
air;  Z.v  is  the  latent  heat  of  vaporization  of  water; 
and  X  is  the  fetch  at  which  the  downwind  measure¬ 
ments  were  made.  Subscript  i  denotes  upwind  (ini¬ 
tial)  conditions,  so  wO  (  and  wq{  are  the  turbulent 
upwind  fluxes.  H^x)  and  H L(x)  are  the  downwind 
surface  fluxes  of  sensible  and  latent  heat. 

Although  the  RPB  upwind  humidity  profiles  con¬ 
sistently  show  a  slight  increase  in  humidity  with 
height  at  the  lowest  levels.  Rider  et  al.  (1965)  ex¬ 
plained  that  this  may  have  been  a  consequence  of  a 
systematic  instrument  error  and  so  reiterated  their 
earlier  conclusion  that  wqj  was  zero.  Then,  also, 
<?j(z)  =  C?j(z')  =  (?j.  Finally,  using  the  integrated, 
two-dimensional  continuity  equation,  we  can  elimi¬ 
nate  W(x,z')  from  eq  13,  which  thus  reduces  to 


r*  — 

/  Hf(x)dx  =  pcpw8<X 


+  Pc, 


■/. 


«h<*> 


U(K2)[T(X,z)-T.{(z)\dz,  (16) 

r  r6q(x) 

|  Hl  (x)  dx  =  pLv  U(X,  z) 

*0  *0 

x  \Q(X,z)-Q{\  dx.  (17) 

Notice  that  the  integrations  now  go  only  to  6h,  the 
thermal  IBL  height,  and  to  5q,  the  moisture  IBL 
height,  since  by  definition  T -  Tt  is  zero  above  6h 
and  Q  -  Q,  is  zero  above  6q. 

Equations  1 6  and  1 7  give  the  total  sensible  and 
latent  heat  exchanged  between  the  grass  surface  and 
the  air  for  the  area  between  x  =  0  and  x  =  X.  I 
obtain  estimates  of  the  average  surface  fluxes  (mW 
cm'2)  by  dividing  the  two  equations  by  the  fetch: 

r6h 

Hs7  =  Wsl  +  PCpX'1  U(T-T)dz,  (18) 
"0 

rq 

HL2=pLvX-'  U(Q  -  Qj)  dz ,  (19) 


,x 

I  |lWs(x)  -pcpwflj)  +  WL(*)}  dx 
■^n 


+  iLv[<?(A',z)-Qi]}  dz 

+  pcp(  [(7(X,z)-(7i(z)l(r.1(z) 
Jo 


-  r,U')]  dz.  (14) 

On  estimating  integrands  in  eq  14,  it  is  evident  that 
the  entire  second  integral  on  the  right-hand  side  is 
negligible  in  comparison  to  the  first  integral  and  I 
henceforth  will  ignore  it. 

Repeating  the  same  analysis  for  the  mass  (water) 
budget  of  the  control  volume  yields 


=  p(  U(X,z)[Q(X,z)-Qj\dz.  (15) 
•'o 


where  Ws1  =pcpw0j. 

After  finding  the  downwind  velocity  profile  using 
the  method  described  in  the  preceding  subsection, 
calculating  the  total  and  average  downwind  heat 
fluxes  from  eq  16-19  is  straightforward.  I  connect 
the  profile  data  with  semilogarithmic  arcs  and  do  the 
indicated  integration  piecewise.  The  only  troublesome 
point  is  assigning  a  surface  roughness  for  the  lower 
limit  of  the  integrations.  I  use  the  value  z02  =  0. 1 4 
cm  given  by  RPB,  setting  f7(z02)  =  0,  T(z02)  =  T02, 
and  (?(zq2)  =  <?02.  The  uncertainties  in  TQ2  and  <?02 
that  I  have  already  discussed  are  not  important  in 
these  integral  calculations.  Because  RPB  measured 
temperature  and  humidity  within  5  cm  of  the  surface, 
errors  in  the  estimated  values  of  T02  or  (?02  could 
lead  to  errors  of  no  more  than  5%  in  the  computed 
fluxes. 

RESULTS  _ 

Energy  budget 

Table  1  fists  the  results  of  my  heat  flux  calculations. 
For  Runs  19  to  43  RPB  included  upwind  and  down¬ 
wind  measurements  of  both  /?n  and  G.  Figure  2 


Table  1 .  Summary  of  the  results  of  the  heat  flux  calculations.  "Up"  refers  to  upwind  of  the 
surface  change,  “down”  to  downwind  of  it.  L  is  the  upwind  Obukhov  length,  q ,  the  downwind 
friction  humidity  (defined  in  eq  38). 


Run 

Fetch 

(m) 

h 

(cm) 

L 

(m) 

u*  (cm/s) 

T  oTC) 

(s/h9) 

Flux  up  (mW  cm~2) 

Flux  down  (mW  cm~^j 

up 

down 

up 

down 

R  n 

G 

H5 

Rn 

G 

Hs 

H/. 

1 

16.0 

_ 

-4.61 

22.2 

28.7 

42.2 

21.6 

-3.06 

56.1 

16.7 

20.4 

67.0 

-14.1 

92.8 

2 

18.5 

- 

-3.86 

21.9 

28.5 

43.9 

21.5 

-3.26 

46.9 

- 

23.4 

59.0 

• 

-11.7 

98.2 

3 

16.0 

- 

-2.68 

19.0 

24.3 

44.2 

21.9 

-3.07 

46.9 

* 

21.8 

59.0 

- 

-14.1 

78.6 

4 

17.4 

320 

-3.92 

21.0 

27.7 

43.4 

30.0 

-3.40 

46.0 

- 

20.4 

52.7 

- 

-28.5 

99.0 

5 

18.1 

- 

-2.66 

17.1 

22.3 

42.4 

22.0 

-2.77 

41.9 

- 

16.2 

52.7 

- 

-20.9 

65.0 

6 

18.S 

238 

-2.39 

15.6 

19.3 

40.1 

19.4 

-2.14 

30.6 

13.5 

41.4 

- 

-8.2 

43.6 

7 

19.5 

221 

-3.41 

16.7 

21.6 

38.8 

23.6 

-1.73 

26.4 

* 

11.7 

33.9 

- 

-19.9 

39.6 

8 

16.5 

166 

-3.17 

17.6 

21.9 

41.1 

25.3 

-1.85 

25.5 

- 

14.9 

33.5 

- 

-13.4 

42.8 

9 

23.5 

273 

-2.01 

13.8 

17.7 

38.5 

23.1 

-1.77 

25.5 

- 

11.3 

33.1 

- 

-4.7 

33.4 

10 

21.9 

218 

-3.58 

16.4 

20.6 

37.6 

23.6 

-1.50 

26.4 

- 

10.6 

33.1 

- 

-4.8 

32.9 

11 

16.6 

233 

-4.67 

17.5 

22.0 

34.5 

22.3 

-1.44 

2.9 

9.9 

8.4 

- 

-22.3 

34.1 

12 

16.8 

166 

-3.36 

15.6 

20.0 

34.9 

21.3 

-1.29 

8.8 

- 

9.8 

15.1 

- 

-15.0 

27.7 

13 

16.0 

189 

-3.86 

18.8 

23.6 

38.4 

26.7 

-2.18 

46.9 

* 

14.7 

50.2 

-16.0 

54.8 

14 

16.2 

172 

-2.84 

17.9 

22.2 

40.8 

24.5 

-2.27 

53.2 

* 

17.4 

61.5 

- 

-10.1 

53.5 

IS 

16.2 

209 

-4.94 

23.2 

29.1 

41.9 

28.4 

-2.09 

50.2 

- 

21.6 

56.9 

- 

-12.0 

64.4 

16 

16.2 

220 

-8.89 

28.2 

36.3 

40.9 

23.7 

-1.97 

44.8 

21.8 

53.6 

- 

-27.2 

75.5 

17 

16.4 

191 

-6.96 

25.0 

32.3 

41.6 

23.1 

-1,51 

41.9 

19.4 

51.5 

- 

-25.6 

51.4 

18 

18.0 

252 

-5.94 

20.8 

26.4 

38.3 

13.8 

-1.52 

16.7 

- 

13.0 

29.3 

- 

-20.0 

42.6 

19 

19.8 

- 

-1.05 

12.9 

15.5 

40.7 

22.2 

-2.40 

33.5 

10.9 

17.6 

43.1 

5.0 

1.0 

39.9 

20 

16.7 

• 

-3.61 

19.1 

24.7 

38.9 

21.5 

-2.84 

35.6 

0.0 

16.4 

44.4 

-2.1 

-14.9 

74.5 

21 

17.9 

- 

-1.92 

14.9 

18.8 

39.6 

26.8 

-2.27 

27.6 

-2.1 

14.7 

31.4 

-2.1 

-5.8 

45.4 

22 

17.8 

- 

-2.43 

17.7 

22.5 

41.3 

27.0 

-2.52 

49.8 

26.8 

19.4 

56.9 

10.0 

1.7 

60.1 

23 

17.9 

- 

-2.19 

17.4 

21.9 

42.2 

23.8 

-2.96 

56.1 

23.0 

20.5 

65.7 

10.9 

-3.1 

69.0 

24 

17.7 

- 

-3.09 

19.4 

24.9 

42.1 

23.1 

-3.46 

57.8 

23.0 

20.2 

67.8 

13.0 

-21.4 

91.2 

25 

17.7 

- 

-2.71 

19.2 

24.7 

43.1 

26.2 

-2.89 

56.1 

23.0 

22.2 

64.9 

12.1 

-0.1 

75.4 

26 

19.1 

- 

-3.80 

20.4 

27.1 

40.6 

19.7 

-2.09 

44.8 

13.8 

19.3 

55.7 

9.2 

-15.1 

60.1 

27 

16.2 

212 

-5.30 

25.7 

34.2 

44.2 

23.3 

-2.01 

41.9 

10.0 

27.3 

53.2 

7.1 

-19.4 

72.2 

28 

17.8 

- 

-2.76 

18.0 

22.9 

41.7 

21.1 

-2.38 

35.6 

10.0 

18.1 

46.5 

7.1 

-6.0 

57.6 

29 

18.1 

236 

-15.17 

30.8 

42.2 

41.1 

26.0 

-2.27 

64.0 

16.7 

16.4 

71.6 

6.7 

-33.7 

99.5 

30 

19.S 

207 

-15.46 

32.6 

42.6 

41.7 

27.3 

-1.60 

64.0 

15.9 

19.2 

71.2 

5.0 

-19.5 

71.3 

31 

16.7 

190 

-14.22 

32.8 

44.0 

42.8 

27.4 

-2.23 

62.0 

15.9 

21.2 

69.9 

5.0 

-37.9 

102.1 

32 

17.8 

210 

-13.78 

31.5 

41.8 

42.3 

28.8 

-2.28 

61.1 

16.7 

19.4 

67.8 

4.6 

-37.9 

99.4 

33 

17.4 

205 

-11.29 

28.3 

37.5 

42.4 

26.9 

-2.33 

59.9 

13.4 

17.2 

67.8 

4.6 

-42.9 

90.9 

34 

17.4 

257 

-8.44 

26.4 

35.3 

43.6 

25.4 

-2.75 

51.9 

10.9 

18.6 

61.5 

4.2 

-36.5 

100.8 

35 

18.3 

191 

-7.73 

24.6 

32.0 

42.6 

22.7 

-2.18 

48.1 

9.6 

16.4 

58.6 

4.2 

-22.9 

72.6 

36 

17.7 

236 

-7.78 

25.9 

33.8 

43,7 

24.9 

-2.72 

46.9 

8.8 

19.1 

56.9 

4.2 

-38.8 

95.2 

37 

22.2 

220 

-7.78 

23.7 

31.1 

42.0 

20.7 

-1.77 

38.1 

7.1 

14.6 

49.4 

4.2 

-18.0 

57.2 

38 

21.5 

194 

-7.06 

21.7 

28.1 

41.4 

19.2 

-1.61 

33.9 

7.1 

12.4 

45.6 

2.9 

-18.0 

47.0 

39 

20.0 

287 

-5.23 

23.0 

28.8 

40.5 

21.1 

-2.34 

44.8 

7.1 

19.8 

54.8 

0.0 

-3.3 

71.3 

40 

16.3 

250 

-7.43 

25.8 

33.7 

39.5 

22.1 

-2.19 

41.0 

7.1 

19.7 

49.8 

0.8 

-15.6 

78.2 

41 

17.7 

228 

-2.60 

17.0 

21.3 

40.6 

21.1 

-2.66 

38.5 

7.1 

16.2 

48.6 

0.0 

-5.0 

59.9 

42 

16.3 

245 

-4.72 

22.4 

28.7 

41.1 

21.1 

-2.40 

35.6 

5.0 

20.2 

46.0 

-0.8 

-18.9 

72.6 

43 

-ILL- 

224 

-4.21 

18.9 

24.4 

38.6 

19.0 

-2.33 

33.5 

5.9 

13.7 

43.5 

0.0 

-10,3 

60.2 
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Figure  2.  Comparison  of  the  computed  fluxes  with  the  measured  values 
of  net  radiation  and  ground  heat  flux. 


compares  these  values  with  the  fluxes  I  compute. 

The  upwind  surface  energy  budget  was  clearly  not 
balanced;  the  downwind  balance  was  better  although 
Hs  +  Hl  seems  somewhat  large.  Taylor  (1971 )  also 
found  the  heat  flux  too  small  to  balance  the  upwind 
budget  for  Runs  7  and  27,  the  only  ones  he  consid¬ 
ered. 

Implicit  in  the  budget  equation  (eq  3)  is  the  steady- 
state  assumption.  The  radiative  surface  temperatures 
RPB  listed  and  the  upwind  surface  temperatures  I 
calculate  suggest  this  may  not  be  a  good  assumption. 
Figure  3  indicates  that  on  26  and  31  January  1961 
the  surface  temperature  rose  to  an  afternoon  maxi¬ 
mum  at  about  1400h  and  then  decreased.  In  other 
words,  the  surface  was  storing  heat  then  releasing  it. 
Approximating  the  storage  term,  we  can  rewrite  eq 
3  as 

Rn-G-Hs-HL-dptct  |f  =  0,  (20) 

where  pe  is  the  density  of  the  tarmac,  ce  is  its  specific 
heat,  and  d  is  the  thickness  of  tarmac  that  undergoes 
heating.  Equation  20  shows  that  the  turbulent  fluxes 
will  be  less  than  Rn-G  when  the  surface  is  heating 
up.  As  the  surface  cools,  this  heat  loss  will  enhance 
the  turbulent  fluxes  and  they  should  exceed  R„  -  G. 


The  upper  traces  in  Figure  3  show  this  energy  con¬ 
version  qualitatively.  Rn-G  -Hs  (WL  was  zero  up¬ 
wind)  has  no  zero  crossing,  although  it  should  if  ex¬ 
cluding  the  storage  term  were  the  sole  reason  for  the 
imbalance.  In  addition,  on  estimating  the  size  of  the 
storage  term  by  estimating  dT/dt  from  Figure  3  and 
using  d  =  1 .5  cm  as  the  thickness  of  the  tarmac  above 
the  heat  flux  plate*,  I  find  that  its  range,  -2  mW 
cm*2,  is  too  small  to  explain  the  15  mW  cm-2  swings 
in  Rn  -  G  -  H%.  There  is  still  an  unexplained  error 
in  the  upwind  heat  budget 

The  assumption  that  //L  was  zero  upwind  cannot 
account  for  this  error.  If  the  humidity  profiles  truly 
increased  slightly  with  height,  as  the  data  show,  they 
would  require  a  downward  flux  of  latent  heat— a  flux 
that  would  increase  the  energy  deficit  rather  than 
decrease  it.  I  feel  that  the  error  is  not  in  the  upwind 
sensible  heat  fluxes  either.  If  these  were  large  enough 
to  balance  the  upwind  budget,  Ws  +  WL  would  be 
substantially  larger  than  Rn  -  G  downwind  (see  eq 
18).  Steady  state  and  the  consequent  balanced  en¬ 
ergy  budget  seems  to  be  a  good  assumption  down¬ 
wind,  however.  The  RPB  radiative  temperatures  do 


•Personal  communication  with  E.F.  Bradley,  CSIRO, 
Canberra,  Australia,  1980. 
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Figure  3.  Temporal  behavior  of  the  calculated  upwind  surface 
temperature  and  of  the  upwind  energy  budget. 


suggest  a  surface  heating  and  cooling  cycle  down¬ 
wind,  especially  on  31  January,  but  no  interrelation¬ 
ship  between  Rn  -  G  -  Hi  -  WL  and  surface  temper¬ 
ature  exists  downwind  as  it  does  upwind.  Although 
the  downwind  surface  temperature  may  have  changed, 
evidently  enough  water  was  available  to  carry  the 
excess  heat  away  without  significant  storage. 

I  finally  must  consider  whether  the  measure¬ 
ments  of  Rn  and  G  could  have  been  in  error. 
Schwerdtfeger  (1976)  discussed  the  difficulties  in 
measuring  ground  heat  flux,  but  the  G  values  tabu¬ 
lated  by  RP8  appear  to  be  the  proper  size  when  com¬ 
pared  with  similar  measurements  summarized  by 
Munn  (1966)  and  Geiger  (1971).  Even  if  die  G 
values  were  systematically  underestimated  by  100%, 
however,  they  would  still  contribute  too  little  to  the 
budget  to  account  for  the  imbalance. 

Because  R„,  on  the  other  hand,  is  the  largest 
term  in  the  upwind  energy  budget,  a  systematic 


error  in  its  measurement  could  well  explain  the  im¬ 
balance.  Using  a  relation  given  by  Seckel  and  Beaudry 
(1973)  to  approximate  total  sun  and  sky  radiation, 
assuming  no  cloud  cover,  and  estimating  incoming 
longwave  radiation  with  Swinbank’s  (1963)  formula 
and  outgoing  longwave  radiation  by  e^oT0]4, 1  can 
estimate  an  upper  bound  on  Rn  over  the  tarmac. 

This  value,  50  mW  cm"2,  is  well  below  some  of  the 
higher  values  RPB  listed  and  thus  seems  to  imply 
that  their  measurements  may  have  overestimated 
the  upwind  value  of  Rn.  Although  correcting  for 
such  an  error  would  bring  the  upwind  energy  budget 
into  better  balance,  I  do  not  feel  that  correction  is 
justified:  estimates  of  incoming  longwave  radiation 
over  short  time  intervals  from  formulas  originally 
intended  for  long-term  averaging  are  too  uncertain 
(Arnfeld  1979).  I  therefore  leave  the  upwind  energy 
budget  as  it  is,  unbalanced,  and  suggest  that  this 
imbalance  may  have  resulted  from  errors  in  the 
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upwind  measurements  of  Rn.  All  other  terms  in 
the  budget  are  of  the  proper  size. 

Latent  heat  flux 

Andreas  et  al.  (1979, 1981 )  measured  sensible 
heat  and  condensate  fluxes  downwind  of  a  step 
change  in  surface  heat  and  moisture  flux.  When  these 
fluxes  are  nondimensionalized  as  Nusselt  numbers 
(Schlichting  1968),  the  Nusselt  numbers  for  both 
sensible  heat  and  condensate  are  equal  for  a  given 
flow  regime.  Because  similar  molecular  processes 
control  the  transfer  of  both  heat  and  moisture 
across  an  interface  and  in  the  light  of  results  reported 
by  Coantic  and  Favre  (1974),  Andreas  et  al.  (1979) 
concluded  that  the  latent  heat  flux  should  also  obey 
the  Nusselt  number  equality.  Figure  4  supports 
their  hypothesis.  In  it  I  have  plotted  the  sensible 
heat  flux  data  of  Andreas  et  al.  (1979),  the  conden¬ 
sate  data  of  Andreas  et  al.  (1981 ),  and  the  downwind 
latent  heat  flux  values  calculated  from  die  RPB 
data.  Although  Andreas  etal.  (1979, 1981)  collected 
their  data  over  arctic  leads  during  winter,  an  environ¬ 
ment  in  almost  total  contrast  to  that  of  the  RPB 
measurements,  one  line  fits  all  three  data  sets;  though 
grossly  dissimilar,  the  different  regimes  obeyed  the 
same  basic  physics.  The  line  in  Figure  4  is 


/V  =  0.08/?x°-76, 

(21) 

where  the  fetch  Reynolds  number. 

Rx  ~  UioqXIv  , 

(22) 

characterizes  the  flow.  The  Nusselt  numbers  for 
sensible  heat  A/$,  latent  heat  /VL,  and  condensate  Nc 
are  defined  as 

Afs  =  //sAf/x(T02-r200) 

=  HfX/pc9D(T02  -  ^20o)> 

(23) 

*L  =  Ht_X/pLvDw (O02  "O200)' 

(24) 

Nc=HcXlpDjC02-C200). 

(25) 

In  eq  22-25  t/200,  ^200> and  C200  are  the  wind 
speed,  temperature,  specific  humidity,  and  conden¬ 
sate  concentration  (in  g/g)  200  cm  above  the  upwind 
surface;  C02  is  the  surface  condensate  concentration 
(Andreas  et  al.  1981);*',  x,  D,  and  Dw  are  the  kine¬ 
matic  viscosity,  thermal  conductivity,  thermal  dif- 
fusivity,  and  water  vapor  diffusivity  in  air;  and  Hc  ; 
is  the  surface  condensate  flux. 


Figure  4.  Nusselt  numbers  from  three  different  data  sets  corre¬ 
lated  with  fetch  Reynolds  number. 
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Figure  5.  Stanton  numbers  plotted  versus  fetch  Reynolds  number. 
The  lines  are  eqs  26  and  27  with  eq  21  substituted  for  N. 


Equation  21  is  slightly  different  than  eq  1,  the 
relation  Andreas  et  al.  (1979)  found  using  only  the 
sensible  heat  flux  data  shown  in  Figure  4.  Numer¬ 
ically,  eqs  1  and  21  are  very  close,  however;  the  two 
are  nearly  coincident  in  the  fetch  Reynolds  number 
region  covered  by  the  sensible  heat  flux  data.  Be¬ 
cause  the  RPB  latent  heat  flux  data  neatly  extend 
this  fetch  Reynolds  number  range,  eq  21  is  the  more 
representative  expression  through  the  entire  region. 

One  might  suppose  that  the  downwind  RPB  sen¬ 
sible  heat  flux  should  also  obey  eq  21 .  It  does  not. 

In  the  RPB  study  the  sensible  heat  flux  behaved  dif¬ 
ferently  than  did  the  fluxes  in  the  three  data  sets 
plotted  in  Figure  4.  Each  of  the  fluxes  represented 
in  Figure  4  went  from  a  near-zero  value  upwind  to  a 
large  positive  value  downwind.  But  during  the  RPB 
experiment  the  sensible  heat  flux  reversed  direction 
within  the  IBL:  it  was  large  and  positive  upwind 
and  large  and  negative  downwind.  The  slope  of  the 
temperature  profile  necessarily  also  reversed  sign 
within  the  IBL.  Under  these  circumstances,  T02  - 
T2 op  would  not  be  the  proper  temperature  scale  in 
eq  23.  There  is  evidently  no  reason,  however,  why 
eq  21  cannot  describe  a  negative  downwind  flux  if 
the  upwind  flux  is  near  zero. 

Stanton  numbers  (Schlichting  1968,  p.  662),  essen¬ 
tially  bulk  aerodynamic  coefficients,  are  sometimes 
used  instead  of  Nusselt  numbers  to  parameterize 
fluxes  in  a  fetch-limited  How  (e.g.  Mangarella  et  al. 
1972, 1973).  Stanton  number:  for  sensible  ( Sts ) 


and  latent  (SfL,  often  called  the  Dalton  number)  heat 
fluxes  are  defined  as 

(26) 

and 

SrL  =  Nl  Kx-’  Sc*1, 

(27) 

where  P  is  the  Prandtl  number,  and  Sc  is  the  Schmidt 
number.  For  a  particular  flow  regime— a  particular 
fetch  Reynolds  number— Ns  =  NL\  hence, 

StL  =  (P/Sc)Sts  =  (DjD)Sts .  (28) 

Because  >  D  (List  1963),  the  Stanton  numbers 
will  not  be  equal.  Figure  5  confirms  this  Stanton 
number  inequality  and  so  emphasizes  the  advantages 
of  using  the  Nusselt  number  as  a  flux  parameter.  The 
lines  in  the  figure  are  eq  26  and  27  with  eq  21  sub¬ 
stituted  for  N.  The  results  reported  by  Mangarella 
et  al.  (1972, 1973)  also  indicate  this  Stanton  number 
inequality.  Although  in  their  results  SfL  is  seldom 
20%  larger  than  Sfs  as  required  by  the  value  of  Dw/D, 
it  is  systematically  larger. 

Surface  stress 

I  have  explained  how  I  found  the  upwind  friction 
velocity  w»up.  After  determining  downwind  velocity 
profiles,  I  computed  the  downwind  friction  velocity 
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Figure  6.  Correlation  between  the  calculated  values  of  up¬ 
wind  and  downwind  friction  velocity. 


4/,  down  by  fitting  a  least-squares  line  through  the 
lowest  three  profile  levels,  U(z 02)  =  0,  U( 27.5),  and 
U( 37.5).  Bradley  (1972)  demonstrated  that  for  pro¬ 
file  levels  near  the  surface,  as  these  are,  the  computed 
value  of  ut  will  be  within  1%  of  the  actual  value  re¬ 
gardless  of  stability. 

Figure  6  compares  downwind  and  upwind  values 
of  u,.  The  result, 

u ‘down  =  1  ,4®  w*up  ~  2.26,  (29) 

which  gives  u*down  in  cm/s  for  u*Up  ,n  cm/s> is  not 
consistent  with  the  RPB  result  w*down/u»up  ~  1.7. 
Equation  29  is  in  good  agreement  with  the  results  of 
Taylor  (1971)  for  Run  7,  atdov/n/u,up  ~  1i43>  and 
of  Raoetal.  (1974)  for  Run  6,  u*do wn/d*up  ~  1.40, 
however.  The  magnitude  of  the  discrepancy  between 
the  RPB  result  and  mine  (i.e.  1.7/1 .4)  is  entirely 
compatible  with  the  difference  in  z0  values  ascribed 
to  the  tarmac.  Previously  I  explained  why  the  zQ 
value  RPB  found  for  the  tarmac,  0.002  cm,  seemed 
too  small.  The  smallness  of  this  value  goes  hand  in 


hand  with  an  underestimate  of  «»up  and  so  with  an 
overestimates  of  w,down/u,up. 

The  implications  of  Bradley’s  (1972)  research  are 
that  we  can  define  a  meaningful  drag  law  regardless 
of  stability  and  even  in  a  developing  flow  if  we  choose 
a  reference  height  low  enough.  Figure  7  shows  that 
drag  law  for  the  downwind  RPB  data, 

w*down  =  0.0646  t/50  +  0.37  .  (30) 

Here  u*down  is  'n  cm/s  w'hcn  ^so< the  wind  *Peed 
50  cm  above  the  surface  at  fetch  X,  is  in  cm/s.  This 
result  is  in  excellent  agreement  with  ones  reported 
by  Andreas  et  al.  (1979)  for  flow  over  arctic  leads, 

u,down  =  0.0643  U50  -  0.77,  (31) 

and  by  Hicks  (1976)  for  more  homogeneous  Aus¬ 
tralian  grassland, 


u,  =  0.0632  U50  +  2.22 . 


(32) 


Figure  7.  Downwind  drag  relation. 


Downwind  humidity  profiles 

One  of  my  principal  reasons  for  undertaking  the 
analysis  of  the  RPB  data  was  an  interest  in  the 
humidity  profile  following  a  step  change  in  surface 
moisture.  We  had  argued  earlier  (Andreas  et  al.  1981 ) 
that  the  condensate  scale  c  in  the  nondimensional 
condensate  gradient  over  arctic  leads, 

0c  =  Wc)~^,  (33) 

i>  a  function  of  two  length  scales,  height  and  fetch. 
When  we  modeled  c  as 

c(X,z)  «  ln[z/A(x)l,  z0<z  <s  ft ,  (34) 

where  h  is  a  fetch-dependent  estimator  of  IBL  height, 
we  found  that  the  condensate  profile  should  have  the 
form 

C(z )K  «  In [z/h) .  (35) 

The  measured  profiles  substantiated  these  heuristic 
arguments  and  also  yielded  the  proportionality  con¬ 
stant  in  eq  35.  We  found 

[-C(e)/c,l*=-0.34lnU//»),  (36) 

with  c*  defined  from  the  surface  value  of  the  con¬ 
densate  flux.  The  excellent  set  of  humidity  profiles 


collected  by  RPB  gave  me  a  chance  to  test  this  model 
on  another  scalar  quantity. 

Figure  8  shows  that  30  of  the  43  RPB  humidity 
profile  sets  follow  an  equation  like  eq  36.  The  line 
in  the  figure  is 


[-AQ/q,}K  =  -0.37  \n(z/h), 

z0<z  <  h, 

(37) 

where  A Q  =  Q(X,z)  -Q.^)  and 

<7*  H^2/pLwk  ^»down  ’ 

(38) 

Hl2  being  the  downwind  latent  heat  flux  calculated 
with  eq  1 9. 

The  multiplicative  constants  in  eqs  36  and  37  are 
slightly  different.  1  suspect  that  the  source  of  this 
difference  is  the  difference  in  z0  values.  The  multi¬ 
plicative  constant  is  evidently  only  a  weak  function 
of  z0/b,  however;  with  the  small  spread  in  the  z0/h 
ratios  of  the  lead  and  RPB  data  sets,  I  can  say  little 
about  its  behavior  except  that  a  conservative  range 
for  it  in  naturally  occurring  internal  boundary  layers 
is  (0.28, 0.42] . 

The  estimate  of  IBL  height  h  was  found  by  plotting 
each  humidity  profile  as  AC? 54  vs  In  /  and  extending 
the  best  straight  line  through  this  profile  to  AQ*4  =  0. 
Figure  9,  which  contains  data  for  both  condensate 
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Figure  8.  Nondimensional  humidity-difference  profiles.  The  30 
profiles  qua!ity<oded  1  or  2  are  shown  here. 


Figure  9.  IBL  height  vs  fetch  for  both  the  conden¬ 
sate  data  of  Andreas  et  al.  (1981)  and  the  30  humid¬ 
ity  profiles  shown  In  Figure  8. 


profiles  (Andreas  et  al.  1981)  and  the  RPB  humidity 
profiles,  shows  that,  as  expected,  h  is  a  function  of 
X.  The  line  is  drawn  by  eye  so  that  h  is  proportional 
toX0-8  (Sutton  1934, 1953,  Calder  1949,  Elliot 
1958,  Philip  1959).  It  is  customary  to  nondimension- 
alize  plots  of  IBL  height  vs  fetch  with  a  surface 
roughness  scale.  The  correlation  in  our  data  does  not 
improve  with  any  such  nondimensionalization,  how¬ 
ever;  the  scaling  lengths  must,  therefore,  be  nearly 
the  same  for  both  data  sets.  Indeed,  the  scalar  rough¬ 
ness  lengths,  zw,  computed  by  the  method  suggested 
by  Garratt  and  Hicks  (1973,  Fig.  1 )  are  0.025  cm 
and  0.023  cm  for  condensate  and  humidity,  respec¬ 
tively.  Choosing  zw  =  0.025,  the  line  in  Figure  9 
becomes 

h/ew  =  1.21  (Xlzw)°*  .  (39) 

The  research  by  Andreas  et  al.  (1981 )  and  this 
study,  are,  to  my  knowledge,  the  only  geophysical- 
scale  experimental  studies  of  IBL  growth  for  passive, 
scalar  contaminants.  Therefore,  although  tne  x-de- 
pendence  of  the  IBL  height  shown  in  eq  39  is  on 


Table  2.  Quality  code*  of  the  AQ 
profiles,  and  wind  direction  during 
each  run.  The  meaning  of  die 
quality  codes  is  given  in  the  text 


Run 

Date 

Quality 

Direction  D 

1 

10  Jan  61 

3 

0 

2 

3 

-30 

3 

3 

0 

4 

1 

-23 

5 

3 

-28 

6 

11  Jan  61 

2 

-30 

7 

1 

~  35 

8 

1 

-14 

9 

1 

-47 

10 

1 

-43 

11 

13  |an  61 

2 

-15 

12 

1 

-18 

13 

2 

0 

14 

1 

-8 

IS 

1 

-8 

16 

1 

-10 

17 

2 

-12 

18 

2 

-27 

19 

26  |an  61 

3 

36 

20 

3 

17 

21 

3 

27 

22 

3 

26 

23 

3 

27 

24 

3 

25 

25 

3 

25 

26 

3 

33 

27 

1 

10 

28 

3 

26 

29 

31  |an  61 

1 

-28 

30 

2 

-35 

31 

2 

-17 

32 

1 

-26 

33 

1 

-23 

34 

1 

-23 

3S 

1 

-29 

36 

1 

-25 

37 

1 

-44 

38 

1 

-42 

39 

11  Feb  61 

1 

-37 

40 

1 

12 

41 

1 

-25 

42 

1 

-12 

43 

1 

22 

sound  footing,  the  multiplicative  constant  must  re¬ 
main  tentative.  Elliott’s  (1958)  model  for  tempera¬ 
ture  as  a  passive  contaminant  and  calculations  by 
Dyer  (1963)-based  on  Philip's  (1959)  theory  of 
local  advection-indicated  a  constant  much  smaller 
than  that  shown  in  eq  39.  In  other  words,  these 


theories  predicted  that  the  IBL  grows  more  slowly 
with  fetch  than  the  experimental  evidence  shows. 

Data  reported  by  Peterson  et  al.  (1976),  however, 
seem  to  support  my  findings;  the  internal  momentum 
boundary  layer  they  studied  grew  higher  than  ex¬ 
pected  on  the  basis  of  numerical  models.  Further 
experimental  work,  with  emphasis  on  a  consistent 
definition  of  IBL  height  and  a  uniform  choice  of 
scaling  length,  is  necessary  to  clarify  these  conflicting 
results. 

Thirty  of  the  43  RPB  humidity-difference  pro¬ 
files  were  proportional  to  (In  z)2.  I  have  quality 
coded  each  profile  with  a  number  from  one  to  three: 

1)  the  profile  is  well  represented  by  AQ^  a  |n  z, 

2)  four  of  the  five  profile  points  show  AQ %  a  In  z, 
and  3)A(?^  «  In  z  is  not  indicated.  Figure  8  con¬ 
tains  the  30  profiles  that  were  assigned  quality  codes 
of  1  or  2.  On  sorting  the  runs  out  as  in  Table  2,  we 
get  some  insight  into  why  all  the  profiles  did  not 
follow  eq  37.  The  13  contrary  profiles  were  collected 
on  only  two  of  the  six  experimental  days,  10  January 
and  26  January.  And  these  profiles  account  for  13 
of  the  15  profiles  collected  on  those  days.  This 
systematic  behavior  suggests  that  an  experimental 

or  calibration  error  on  each  of  these  two  days  could 
have  masked  the  (In  z)2  behavior  of  the  humidity- 
difference  profiles.  If  only  two  of  the  20  thermo¬ 
couples  RP8  used  to  measure  upwind  and  downwind 
wet  and  dry  bulb  temperatures  had  bias  errors,  the 
AQft  «  In  z  profile  form  would  have  been  lost.  The 
evidence  that  such  errors  occurred  on,  at  most,  two 
days  is  testimony  to  the  experimental  care  Rider, 
Philip,  and  Bradley  took. 

Table  2  also  suggests  that  the  quality  of  the  humid¬ 
ity  profiles  may  have  been  influenced  by  wind  direc¬ 
tion.  Or.  26  January— when  nine  of  the  ten  profiles 
were  quality  3-the  wind  direction  was  positive  (0° 
is  along  the  site  axis).  During  only  two  of  the  re¬ 
maining  33  runs  was  the  wind  direction  positive. 

This  is  strong  circumstantial  evidence  that  the  flow 
was  not  in  equilibrium  with  the  tarmac  when  its  dir¬ 
ection  was  ~  30s  but  was  apparently  disturbed  by 
the  large  buildings  upwind  in  that  direction.  RPB 
explained,  in  fact,  that  even  when  the  wind  direction 
was  nearer  0°,  these  buildings  were  the  source  of 
"odd  gusts."  I  feel,  therefore,  that  flow  distortion 
by  these  buildings  voided  the  upwind-equilibrium 
assumption  on  26  January,  resulting  in  the  poor 
quality  of  the  humidity  profiles  on  that  day. 


DISCUSSION 

Results  in  the  last  section  present  a  paradox.  The 
Nusselt  numbers  for  three  different  scalar  fluxes, 


sensible  heat,  latent  heat,  and  condensate,  obey  the 
same  fetch  Reynolds  number  relation,  but  tempera¬ 
ture  profiles  over  arctic  leads  do  not  have  the  same 
form  as  the  condensate  and  RPB  humidity  profiles. 
One  set  of  results  suggests  that  scalar  quantities  are 
transferred  similarily  in  an  IBL,  while  a  second  set 
of  results  shows  inherent  differences  in  the  transfer. 
On  the  basis  of  current  boundary  layer  theory,  these 
results  are  not  really  incompatible,  however.  The 
scalar  quantities  we  have  studied  arc  transferred 
across  a  surface  by  virtually  identical  molecular 
processes-which  transfer  determines  the  Nusselt 
number.  But  they  need  not  obey  similar  physics 
in  the  turbulent  IBL,  where  the  profiles  were  meas¬ 
ured. 

Although  I  have  computed  the  downwind  surface 
fluxes  from  the  measured  profiles  according  to  eq  1 8 
and  19,  eq  16, 17  and  23-25  show  that  the  surface 
flux  of  a  scalar  5  is  related  to  the  Nusselt  number  by 

/Vj  =  (S0  -  ■S200)*1  dx . 

(40) 

Here  Os  is  the  molecular  diffusivity  of  the  scalar; 

50,  its  surface  value;  5200,  its  value  200  cm  above 
the  upwind  surface;  and  Hs  is  now  the  kinematic  sur¬ 
face  flux  (i.e.  with  constants  such  as  pcp  and  piv 
omitted). 

Hs  can  be  computed  from  the  S  gradient  through 
the  molecular  sublayer, 


profiles  of  the  other  scalars  are  also  only  weak  func¬ 
tions  of  their  molecular  diffusi vi ties — as  is  tempera- 
ture-and  if  the  scaling  performed  in  the  preceding 
paragraph  is  appropriate  for  all  the  scalars,  eq  42 
implies  that  the  corresponding  Nusselt  numbers  will 
be  nearly  identical  functions  of  the  fetch  Reynolds 
number  in  the  IBL.  The  Nusselt  number  equality 
that  we  have  found  is  a  posteriori  evidence  for  the 
validity  of  these  assumptions;  there  is  similarity  in 
the  transport  of  scalars  through  a  surface  molecular 
sublayer.  The  predicted  dependence  of  N  on  molec¬ 
ular  diffusivity  is  such  a  small  effect  that  it  is  lost 
amid  the  experimental  uncertainties. 

In  the  turbulent  IBL,  on  the  other  hand,  the 
transfer  processes  need  not  be  the  same  for  each 
scalar.  Warhaft  (1976)  has  developed  a  simple  model 
for  turbulent  transfer  in  a  horizontally  homogeneous 
surface  layer.  Although  his  model  is  not  directly 
applicable  to  internal  boundary  layers,  it  does  indi¬ 
cate  the  likelihood  of  differences  in  the  turbulent 
transfer  of  scalars:  differences  that  are  here  our 
concern.  We  can  thus  gain  some  insight  from  his 
model.  The  equations  that  Warhaft  (1976)  derived 
for  the  turbulent  fluxes  of  temperature  and  humidity 
are 

^0  =  -  'A{glw2)(o2IT0  +  0.61  0^j| 

(43a) 

and 


Hs(x,z-0)  =  -Ds 


lsi*A 

d z 


(41) 


Schlichting  (1968,  p.  284)  showed  that  a  scaling 
length  for  7  in  fetch-limited  boundary  layers  is 
(vx/U 200)^  •  Such  a  scaling  length  seems  appropriate 
in  Nusselt  number  computations  for  any  of  the 
scalars  we  are  discussing,  because  for  large  Rx  the 
Nusselt  number  cannot  depend  on  stability  (Andreas 
1980).  With  this /-scale  in  eq  41,  eq  40  becomes 


a5  (x,/,/?x,wps) 
£ 


(42) 

where  the  wavy  overbars  indicate  new  nondimensional 
quantities:  x  =  x/X,  z  =  (/2£/2oo/Wf)>*> an<*  ^ 

=  (S  -  5  200)  /  (Sq  •  Sjqq). 

Schlichting  (1968)  went  on  to  show  that  in  an 
IBL  the  nondimensional  temperature  profile  evalu¬ 
ated  at  the  surface,  df/bz  t^»0,  is  a  weak  function 
of  Prandtl  number  (~/,,  /3).  If  the  nondimensional 


iW7=  ~K\&  '  *(^2K°-61  q2  +  */r0j|  . 

(43b) 

Here  T0  is  a  temperature  representative  of  the  sur¬ 
face  layer;  and  K  =  Vi  (w2  e/eA),  where  F  (=£ 

+  v?+  ^  is  the  turbulent  kinetic  energy,  €  is  the 
energy  dissipation  rate,  and  A  is  a  dimensionless 
constant 

I  nondimensionalize  eq  43  with  the  standard  de¬ 
viations  of  temperature,  og  -  6^Vl,  humidity,  oq 
=  q2V> , and  vertical  velocity,  o*.  =  w2  Vll  and  with 
the  temperature-humidity  correlation  coefficient, 
pgq=0qlo9oq: 

*tilo*og  =  -(gKloJ)\{ov/2lgoe)  ~ 

-  'M°0IT0  +  0.61  oq  p9p)l  ,  (44a) 

*V/°w°q  =  -  (^AC/cr^3 )  [  (aw2/sr<zq )  ~ 

-Moe  p9q/r0  +  0.61  oq)l  . 


(44b) 


If  p0q  =  I,  these  two  equations  have  identical  non- 
dimensional  forms.  The  fluxes  then  clearly  result 
from  two  sources,  a  flux  down  the  scalar  gradient 
and  a  nongradient  component  associated  with  the 
turbulence.  Notice  that,  depending  on  its  sign,  the 
scalar  gradient  can  drive  an  upward  or  a  downward 
flux;  but  the  nongradient  flux  is  always  upward- 
hence,  at  times  it  is  countergradient  (Deardorff  1972). 

On  rearrangement,  eq  44a  and  44b  yield  the  non- 
dimensional  gradients: 

-(<’w2/gag}~=  crw2  »«lgKoe  -  'Moe/T0  +  0.61  oq  p0q)  , 

(la)  (2a)  (3a)  (45a) 

0  .30  o  — 

-(»*  /8° q>^7  =  <V  wqtgKoq  -Vi|o0  PeqITQ  +  0.61  oq)  . 

Ob)  (2b)  (3b)  (45b) 

From  these  it  is  evident  that  the  scalar  profiles  are 
determined  by  both  the  boundary  conditions  (terms 
2a  and  2b)  and  the  turbulence  within  the  layer  (terms 
3a  and  3b).  The  nondimensional  profiles  will,  there¬ 
fore,  be  dissimilar  if  any  one  of  the  following  is  true: 

1 .  If  the  boundary  conditions  are  different— for 
example,  if  wq  is  positive  and  wO  is  negative. 

2.  If  the  temperature  and  humidity  fluctuations 
are  large  but  poorly  correlated. 

3.  If  the  buoyancy  effects  induced  by  the  tur¬ 
bulence  are  negligible  in  either  eq  45>  or  45b, 
i.e.,  if  either  term  3a  or  3b  can  be  ignored,  but 
not  both. 

Experimental  work  by  McBean  and  Miyake  (1972), 
McBean  (1973),  and  Friehe  et  al.  (1975)  corroborates 
these  theoretical  conclusions. 

The  humidity  variance  is  rarely  large  enough  to 
contribute  to  the  buoyancy  term'  (3a  or  3b)  in  eq 
45a  and  46b  (Friehe  et  al.  197  humidity  is  usually 
a  passive  scalar  with  respect  to  the  transports  of  both 
temperature  and  humidity.  Temperature  is  also 
passive  with  regard  to  humidity  if  (o0q  I  is  small  and 
oe  is  not  extremely  large.  Only  when  ogl2T0  is  an 
appreciable  fraction  of  term  la  does  temperature 
act  as  an  active  scalar  in  the  temperature  equation. 

These  criteria  for  judging  whether  scalars  are  active 
or  passive  extend  those  usually  seen  in  the  literature. 

For  example,  McBean  and  Miyake  (1972)  wrote,  “A 
passive  scalar  is  a  scalar  whose  variations  do  not  sig¬ 
nificantly  affect  the  buoyancy  of  a  parcel  of  air.  An 
active  scalar  is  a  scalar  whose  variations  do  affect  the 
buoyancy."  Through  eq  45a  and  45b  I  have  shown, 
however,  that  a  particular  scalar  can  contribute  to 
the  buoyancy  transfer  of  one  scalar  quantity  yet  not 
to  another.  Such  a  scalar  is  neither  totally  active  nor 
totally  passive. 
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Although  we  must  be  cautious  in  drawing  infer¬ 
ences  about  scalar  transfer  in  internal  boundary  layers 
from  eq  45a  and  45b,  their  buoyancy  and  gradient 
terms,  at  least,  cannot  be  wrong.  The  equations  thus 
provide  some  insight  into  the  profile  results  discussed 
in  the  Downwind  humidity  profiles  section.  The  RPB 
temperature  and  humidity  profiles  within  the  IBL 
will  clearly  be  different  because  the  boundary  condi¬ 
tions  are  different;  wd  is  downward  at  the  grass  sur¬ 
face  and  upward  near  the  top  of  the  IBL,  while  wq 
is  upward  everywhere.  p0q  will  change  from  negative 
to  positive  as  the  height  increases,  and  I  suspect  its 
absolute  value  will  consequently  be  small.  Term  3b 
should  thus  be  negligible-buoyancy  does  not  affect 
humidity  transfer  for  the  RPB  data. 

Condensate  over  arctic  leads  should  obey  a  trans¬ 
fer  equation  just  like  that  for  humidity.  Andreas  et 
al.  (1981)  have  shown  that,  because  of  the  necessity 
of  diffusion  through  a  molecular  boundary  layer 
surrounding  a  condensate  droplet,  droplets  over  leads 
change  little  in  size  for  several  seconds  despite  step 
changes  in  ambient  humidity  or  temperature.  Over 
leads,  away  from  the  immediate  vicinity  of  the  sur¬ 
face,  condensate  is  thus  a  quasi-conservative  property 
and  will  obey  the  same  diffusion  equation  as  tempera¬ 
ture  and  humidity  (Pasquill  1974).  A  condensate 
transfer  equation  analogous  to  that  for  humidity  (eq 
45b)  then  follows  directly.  But  again  the  buoyancy 
term  will  likely  be  negligible.  Although  the  tempera¬ 
ture  variance  over  leads  can  be  extreme -0.6  to  1 ,0°C2~ 
the  correlation  between  temperature  and  condensate 
should  be  small,  since  a  positive  temperature  fluctua¬ 
tion,  for  example,  would  lower  the  probability  of 
formation  of  a  condensate  droplet  by  lowering  the 
relative  humidity.  And  I  have  already  mentioned  how 
slowly  the  small  condensate  droplets  respond  to 
changes  in  humidity  and  temperature.  Consequently, 
as  with  humidity  during  the  RPB  experiment,  buoy¬ 
ancy  does  not  affect  condensate  transfer:  both  scalars 
should  obey  the  same  physics. 

For  temperature  and  humidity  over  leads,  in  con¬ 
trast,  the  buoyancy  term  contributes  several  percent 
to  both  transfer  equations  because  of  the  large  value 
of  o9  and  the  high  temperature-humidity  correlation 
over  the  warm,  evaporating  surface  (e.g.  Phelps  and 
Pond  1971,  Friehe  etal.  1975;Wesely  and  Hicks 
1978).  Thus,  though  the  boundary  conditions  on  w6, 
wq,  and  wc  over  leads  and  on  wq  in  the  RPB  experi¬ 
ment  are  the  same,  the  corresponding  scalar  profiles 
will  not  all  be  similar.  Temperature  and  humidity 
profiles  should  scale  similarly,  and  the  condensate 
and  RPB  humidity  profiles  should  scale  similarly, 
but  the  two  groups  must  differ  because  the  latter  are 
unaffected  by  buoyancy. 


CONCLUSIONS 


The  well-respected  data  set  of  Rider  et  al.  (1963) 
is  a  standard  with  which  to  test  models  of  flow 
following  a  step  change  in  surface  conditions.  For 
each  run  in  the  data  set  I  have  calculated  several 
quantities  that  have  not  been  reported  before.  Prin¬ 
cipal  among  these  are  the  sensible  and  latent  heat 
fluxes  and  friction  velocity  upwind  and  downwind. 

The  upwind  surface  energy  budget  is  not  balanced 
in  the  steady-state  form.  Including  a  storage  term 
qualitatively  improves  the  balance,  but  there  is  still 
a  discrepancy  that  may  be  due  to  an  error  in  the  up¬ 
wind  measurement  of  net  radiation.  The  downwind 
energy  budget  is  in  approximate  steady-state  balance. 

The  surface  latent  heat  flux,  which  goes  from 
virtually  zero  upwind  to  a  large  upward  value  down¬ 
wind,  is  related  to  the  fetch  Reynolds  number.  When 
the  latent  heat  is  nondimensionalized  as  a  Nusselt 
number, 

A'l  =0.08/?x°-76  .  (46) 

This  expression  is  essentially  the  same  as  one  found 
to  describe  sensible  heat  and  condensate  fluxes  over 
arctic  leads.  I  conclude  that  all  Nusselt  numbers  non- 
dimcnsionalizing  scalar  surface  fluxes  following  a 
step  increase  in  the  flux  are  equal  for  a  given  flow 
regime. 

The  RPB  humidity-difference  profiles  are  well- 
represented  by 

{-AQIq*\'A  =-0.37  \n(zjh),  (47) 

where  h  is  a  fetch-dependent  indicator  of  internal 
boundary  layer  height.  This  form  of  the  humidity- 
difference  profile,  with  the  exception  of  the  multi¬ 
plicative  constant,  is  identical  to  that  found  for  con¬ 
densate  profiles  at  the  downwind  edge  of  arctic  leads. 
The  common  links  between  these  two  scalars  are 
that  they  have  similar  boundary  conditions  and 
neither  is  affected  by  buoyancy  transfer.  I  recommend 
looking  for  this  relationship  in  the  profiles  of  other 
such  scalars  in  fetch-limited  flows. 

Equations  46  and  47  imply  two  distinct  methods 
of  estimating  the  latent  heat  flux  in  fetch-limited 
flows  from  easily  measured  quantities.  Equation  46 
calls  for  the  following  five  measurements-  surface 
temperature  (for  computing  Q02,  p,  v,  Lv,  and  Dw), 
wind  speed  at  2  m,  upwind  humidity,  fetch,  and 
atmospheric  pressure  (also  needed  to  find  p).  Seven 
measurements  are  necessary  for  determining  the  flux 
from  eq  47:  upwind  humidity  at  two  heights,  down¬ 
wind  humidity  at  the  same  two  heights,  surface 


temperature  (for  computing  p  and  Lv),  wind  speed 
50  cm  above  the  downwind  surface  (which  gives  u, 

[ eq  30) ),  and  pressure  (again  for  p).  Although  this 
second  method  requires  two  more  measurements 
than  the  first,  it  yields  two  more  quantities,  u,  and 
the  IBL  height  h. 

Some  of  my  results  are  not  consistent  with  the 
data  reported  in  RPB  and  so  raise  questions  about 
the  accuracy  of  their  data.  The  most  striking  example 
of  possible  inaccuracies  is  in  the  humidity  profiles 
of  1 0  january  and  26  January .  As  I  have  explained, 
there  was  evidently  either  a  calibration  error  on  those 
days  or  the  flow  was  disturbed  by  buildings  upwind. 
The  upwind  value  of  z0  (0.002  cm)  RPB  gave  is  very 
small.  I  feel  it  is  too  small;  using  it  I  could  not  re¬ 
trieve  reasonable  upwind  and  downwind  velocity  pro¬ 
files  from  the  average  profiles  RPB  listed.  A  value 
of  z0  =  0.01  cm  upwind  leads  to  upwind  and  down¬ 
wind  profiles  that  compare  well  with  the  few  indi¬ 
vidual  profiles  available  and  is  more  compatible  with 
aerodynamic  considerations.  The  value  of  1 .7  RPB 
quoted  for  the  ratio  u,down/u»up  is  significantly 
larger  than  the  value  of  1 .4  I  would  give  it  on  the 
basis  of  eq  29.  This  disparity  is  a  consequence  pri¬ 
marily  of  the  larger  zQ  value  I  ascribe  to  the  upwind 
tarmac.  Lastly,  I  compute  downwind  surface  temper¬ 
atures  lower  than  RPB  measured.  My  values  are  con¬ 
sistent  with  the  downwind  temperature  and  humidity 
profiles,  while  the  surface  temperatures  of  RPB  imply 
a  complex  shape  for  the  temperature  profile  and  a 
grass  surface  well  below  saturation. 
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